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LONG-TERM  GOALS 

The  overall  long-term  goal  for  this  project  is  to  determine  if  and  how  acoustic  coherent  backscatter 
enhancement  (ACBE)  can  be  used  for  classification  of  active  sonar  returns  in  a  wide  variety  of  ocean 
environments.  During  its  first  three  years,  this  project  has  focused  on  simulations  of  acoustic  multiple 
scattering  from  two-  and  now  three-dimensional  aggregations  of  omni-directional  point  scatterers  to 
determine  the  parametric  realms  in  which  ACBE  might  be  observed,  and  its  characteristics  when  it  is 
observed. 

OBJECTIVES 

The  detailed  objectives  of  the  current  research  effort  are  to  determine  the  parametric  dependence  of 
ACBE  peak  amplitude,  peak  emergence  rate  as  the  number  of  observations  increases,  peak  angular 
width,  and  peak  time  dependence.  Here  the  independent  parameters  are  the  range  between  the 
scattering  aggregation  and  the  receiving  array,  the  receiving  array  characteristics,  incident  wave 
characteristics  (wave  front  shape,  waveform,  frequency,  bandwidth,  duration),  and  aggregation 
characteristics  (scatterer  cross  section  and  mean  spacing,  overall  aggregation  size  and  shape,  etc.). 
Eventually,  underwater  waveguide  characteristics  will  be  considered  as  well. 

APPROACH 

The  current  approach  involves  numerical  evaluation  of  the  fundamental  equations  of  multiple 
scattering  from  an  aggregation  of  omni-directional  point  scatterers1.  If  y/{ r)  is  the  harmonic  acoustic 
pressure  field  at  frequency  co  at  the  point  r  and  y/o(r)  is  the  harmonic  field  incident  on  the  aggregation 
of  scatterers  located  at  r„,  then 

N 
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where  if/ [pc)  is  the  scattered  field  and  is  given  by  the  sum  in  (1),  A  is  the  number  of  scatterers,  gn  is  the 
scattering  coefficient  of  the  nth  scatterer,  y/n( r„)  is  the  field  incident  on  the  nth  scatterer, 

wM=y/0(O+  Z  gjVj(rj)G{*nJj)  (2) 
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and  G(rn,  ry)  is  the  free-space  Green's  function  between  the  locations  r„  and  r y, 


G(r„,ry)  = 


exp{-^0r„-r/} 


where  &0  =  cole  is  the  wave  number  magnitude  of  the  incident  field,  c  is  the  sound  speed,  and  i  = 
When  the  incident  field  and  the  scattering  coefficients  are  known,  (2)  can  be  written  N  times,  once  for 
each  scatterer  1  <  n  <  N,  and  these  N  algebraic  equations  can  be  solved  to  determine  y/n(r „).  The  total 
field  at  any  location  is  then  recovered  from  (1)  and  (3)  using  the  known  yfo(r),  the  known  g„,  and  the 
calculated  y/n(r „).  This  formulation  is  akin  to  the  direct  boundary-integral  formulation  of  computational 
acoustics  with  one  computational  element  assigned  to  each  scatterer.  The  computational  burden  of  this 
approach  is  set  by  the  inversion  of  the  fully-populated  N-by-N  algebraic  system  that  determines 


In  the  current  investigation,  the  scatterers  are  placed  randomly  with  an  average  spacing  s,  and  are 
identical  without  internal  lossless  so  the  N  scattering  coefficients  are  all  the  same 
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where  <js  is  the  scattering  cross  section.  In  the  current  investigations,  as  is  considered  an  independent 
parameter  within  the  constraint  imposed  by  conservation  of  acoustic  energy:  os  <  An  fc2  =  A2  jn.  For 
the  current  calculations,  the  incident  field  y/, o(r)  is  a  plane  wave  with  amplitude  A  and  wave  number 
vector  k0,  and  the  backscatter  direction  (</>  =  cp  =  0)  is  defined  with  respect  to  k0  as  shown  in  Figure  1. 

Three  dimensional  rectangular  aggregations  of  scatterers  with  depth  X,  width  Y,  and  height  Z  have 
been  the  primary  geometry  considered  this  fiscal  year  because  this  geometry  readily  allows  both  near- 
field  and  far-field  investigations  with  a  linear  receiving  array  of  length  L  when  the  number  of  scatterers 
is  limited  to  a  few  thousand. 


To  search  for  the  presence  or  absence  of  ACBE,  the  scattered  field  y/s(r)  predicted  by  (1)  is  calculated 
at  the  elements  of  the  receiving  array,  and  aggregation-array  distance  R  is  varied  to  put  the  receiving 
array  in  the  near-field,  where  the  array's  beam-steering  angle  (j)  is  relevant,  or  in  the  far-field  of  the 
aggregation,  where  the  azimuthal  scattering  angle  (p  is  relevant.  For  near-field  calculations,  the  results 
are  presented  as  B( (/))l[B]exciluuno peak  vs.  (/),  where  B{(j>)  is  the  beamformed  output  of  the  receiving  array, 
and  [B]exciudingPeak  is  the  array's  average  beamformed  output  in  directions  near  backscatter  but  excluding 
the  ACBE  peak.  The  independent  parameters  of  these  investigations  are  A,  k(h  crs,  s,  R,  L,  or  cp ,  X,  Y, 
and  Z,  but  the  current  problem  can  be  stated  in  dimensionless  terms  by  normalizing  field  values  with 
the  incident  wave  amplitude  A,  and  by  using  wave  number-scaled  lengths.  This  reduces  the  number  of 
independent  parameters  on  which  ACBE  may  depend  from  ten  to  eight. 

The  current  focus  of  this  research  effort  is  on  determining  and  understanding  the  parametric 
dependence  of  the  ACBE  peak's  height  and  width  in  the  near  field  on  the  eight  dimensionless 

1  /O 

parameters  (ko<rs  ,  kos,  koR,  koL,  4>or  cp,  leX,  koY,  and  koZ).  In  particular,  the  current  simulations  show 
the  ACBE  peak  width  is  set  by  the  dimensionless  length  koL  of  the  receiving  array,  while  the  peak 
enhancement  depends  on  a  combination  of  /i0crvl/2,  k0s,  and  the  beam  steering  angle  between  the  center 
of  the  receiving  array  and  the  edge  of  the  rectangular  aggregation,  <f>edge  [=  tan  1  f  T72Z2)] . 


2 


These  ACBE  investigations  are  the  current  doctoral  research  of  Ms.  Adaleena  Mookerjee.  She  is  a  US 
Citizen  and  a  Ph.D.  candidate. 


Figure  1.  Geometrical  configuration  of  the  multiple  scattering  simulations.  A  plane  wave  with 
amplitude  A  and  wave  number  vector  ka  impinges  on  a  rectangular  aggregation  of  scatterers  with 

dimensions  X,  Y,  and  Z.  Here  (p  is  the  scattering  angle  defined  from  the  origin  of  coordinates,  and  (j) 
is  the  beam-steering  angle  defined  from  the  broadside  direction  of  the  receiving  array  of  length  L 
that  lies  a  distance  R  from  the  aggregation.  The  backscatter  direction  is  cp=  <f>=  0. 

WORK  COMPLETED 

In  past  years,  the  simulation  code  has  been  shown  to  properly  recover  Bragg  scattering  angles  from 
uniformly  spaced  scatterers,  and  to  conserve  acoustic  energy  to  better  than  ±1%  of  the  acoustic  power 
incident  on  one  scatterer  even  when  thousands  of  scatterers  are  present.  To  complete  the  verification 
and  validation  of  the  simulations  two  more  tests  were  undertaken:  (/')  matching  of  the  simulated 
effective  wavenumber  within  the  aggregation  to  Foldy's  effective  medium  theory,  and  (if)  showing 
consistency  of  the  simulated  ACBE  peak  heights  and  widths  with  optics  experiments.  These  tests  were 
completed  and  an  empirical  coherence  function  for  CBE  was  found  along  the  way.  In  addition,  the 
parametric  dependence  of  the  near-field  CBE  peak  height  on  ko<7sm  and  kos  is  now  under  investigation. 
The  results  are  summarized  in  the  next  section  and  in  a  soon-to-be-submitted  manuscript2. 

RESULTS 

Sample  effective  medium  wave  number  matching  results  are  shown  in  Figure  2.  Here  the  real  part  of 
the  ensemble-average  acoustic  pressure  (black  curve)  inside  a  random  aggregation  of  point  scatterers  is 
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plotted  along  with  the  real  part  of  the  incident  harmonic  plane  pressure  wave  (red  curve)  as  a  function 
of  the  incident-wave-number-scaled  distance,  kyc,  inside  the  aggregation.  Here,  the  simulated  wave 
number  (wavelength)  is  slightly  higher  (lower),  and  this  wave  number  change  is  predicted  by  Foldy's 
effective  medium  theory;  kc2  =  k02  +  4^Vv|gi|,  where  ke  =  2 nlXe  is  the  wave  number  inside  the  scattering 
aggregation,  Ae  is  the  wavelength  inside  the  scattering  aggregation,  Nv  is  the  number  of  scatterers  per 
unit  volume,  and  g\  is  given  by  (4).  For  the  conditions  of  these  simulations,  the  ratio  kjk  is  predicted  to 
be  1.0182,  and  the  value  for  this  ratio  determined  from  the  simulations,  1.0183,  is  in  good  agreement 
with  the  value  from  theory. 


Figure  2.  A  comparison  of  the  real  part  of  the  ensemble-average  harmonic  pressure  (black  curve) 
inside  a  random  aggregation  of point  scatterers  and  the  real  part  of  the  incident  plane  wave  (red 
curve).  The  horizontal  axis  is  the  distance  scaled  by  the  wave  number  of  the  incident  plane  wave. 
The  wavelength  inside  the  aggregation  is  slightly  shorter,  and  simulation  results  shown  here 
conform  to  the  theoretically-expected  wavelength  change. 

As  a  final  validation  test,  results  from  the  simulations  were  compared  to  optical  CBE  experiments  . 
This  proved  challenging  for  several  reasons.  The  optics  experiments  involved  polarized  green  laser 
light,  trillions  of  0.1 09-pm-diameter  plastic  beads,  and  an  angular  resolution  of  0.25  milli-radians 
(mrd)  while  the  simulations  involved  longitudinal  acoustic  waves,  3900  point  scatterers,  and  an 
angular  resolution  of  ~2°  (30-40  mrd).  To  show  that  the  results  from  the  experiments  and  simulations 

1  /O 

were  comparable,  scattering  cross  sections  (kocrs  =  0.23,  0.21),  average  scatterer  spacings  (kos  =  3.1, 
3.2),  and  mean-free  paths  ( kf  =  925,  926)  were  closely  matched.  To  account  for  the  difference  in 
angular  resolution,  a  single  analytical  coherence  function  was  identified  and  used  to  determine  the 
shape  of  the  CBE  peak  using  beamforming4  when  the  receiving  array  length  is  adjusted  to  alternately 
match  the  resolution  of  the  optical  experiments  and  the  simulations.  This  coherence  function  is: 
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where  Ay  is  the  spatial  shift  variable,  and  Ad,  A /,v,  and  <ry  where  chosen  to  produce  the  results  shown  in 
Figure  3.  The  green  curves  in  both  panels  of  Figure  3  come  from  the  coherence  function  and  both 
appropriately  match  the  corresponding  experimental  and  simulation  results.  Thus,  the  simulations  and 
the  experiments  -  although  different  in  resolution  -  are  consistent  with  each  other. 
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Figure  3.  Ensemble  averaged  beamformed  intensity  ratio,  B( </j)/B ( <fi)  excluding peak,  from  the  optics 
experiments  (left)  and  from  the  simulations  (right)  versus  angle  (f),  where  <f>=  0  is  the  exact 
backscatter  direction.  The  important  difference  between  the  two  panels  is  the  difference  in  the 
angular  resolution  of  the  receiving  array  (0.25  mrd  on  the  left  and  2°  on  the  right).  In  both  panels, 
the  green  curve  is  obtained  from  (5)  using  the  same  values  of  Ad,  A/,s,  and  crr  Thus,  the  results  are 
consistent  even  though  the  angular  resolutions  differ  by  approximately  two  orders  of  magnitude. 


Overall,  after  four  verification  and  validation  checks,  no  deficiencies  have  been  found  in  the 
simulations.  However,  while  reaching  the  results  shown  in  Figure  3,  the  reason  for  the  prior  findings  of 
much  larger  than  expected  enhancements  has  been  traced  to  the  finite  size  and  nominally  flat  front  of 
the  rectangular  aggregation  of  scatterers  used  in  the  simulations  and  depicted  in  Figure  1 .  Current  work 
involves  production  of  a  functional  scaling  for  the  peak  ratio  at  <f>=  0  in  the  near  field  of  the 
aggregation  based  on  simulation  results  for  different  combinations  of  ko<jsm  and  kos. 

IMPACT/APPLICATION 

In  broad  terms,  this  project  ultimately  seeks  to  determine  if  and  how  ACBE  might  be  exploited  for 
active  sonar  applications.  In  particular,  if  successful,  it  should  prove  useful  for  remote  classification, 
because  a  large  sonar  return  from  a  single  large  scatterer  will  likely  not  display  ACBE  while  a 
similarly  large  sonar  return  from  an  aggregation  of  many  small  scatterers  may  display  ACBE.  Thus, 
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this  research  effort  may  eventually  impact  how  active  sonar  signals  are  processed  and  displayed  for 
tactical  decision-making  related  to  classification. 

TRANSITIONS 

The  results  of  this  research  effort  should  aid  in  the  design  of  active  sonar  signal  processors  for  tactical 
decision  aids.  However,  at  this  time  no  direct  transition  links  have  been  established  with  more  applied 
research  programs.  Once  the  current  simulation  capability  has  been  more  fully  exploited  for  parameter 
studies,  and  promising  results  have  been  obtained,  a  transition  path  through  NRL  or  one  of  the  Navy's 
Warfare  Centers  will  be  sought. 

RELATED  PROJECTS 

This  project  is  related  to  the  other  projects  funded  under  ONR’s  2010  basic  research  challenge 
program.  In  particular,  the  efforts  by  Prof.  Feuillade  in  Chile  and  Prof.  Sabra  at  Georgia  Tech  are  most 
closely  related. 
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